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Abstract: We derive general tree- level recursion relations for amplitudes which include 
massive propagating particles. As an illustration, we apply these recursion relations to 
scattering amplitudes of gluons coupled to massive scalars. We provide new results for 
all amplitudes with a pair of scalars and n < 4 gluons. These amplitudes can be used as 
building blocks in the computation of one-loop 6-gluon amplitudes using unitarity based 
methods. 
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1. Introduction 



The twistor string description of N = 4 Yang-Mills proposed by Witten [1] has been 
an inspiration for much of the recent progress in detailed calculations of multi-particle 
amplitudes in gauge theory. New ideas have led to the development of powerful new 
formalisms for these calculations, most notably the MHV rules of Ref. [2], the Britto- 
Cachazo-Feng (BCF) recursion relations of Refs. [3, 4], and the generalized unitarity cuts 
in the complexified Minkowski space of Ref. [5]. Applications of these new formalisms 
together with the classic unitarity based approach of Refs. [6, 7] have led to a dramatic 
progress in calculations of amplitudes. 

At tree level, new and compact results for scattering amplitudes were derived using the 
MHV rules [2, 8-18], and the recursion relations [3, 4, 19-23]. External massive particles, 
such as Higgs bosons [14, 16] and electroweak bosons [15], have also been included using 
generalizations of MHV rules. On the other hand, the BCF recursion relations have been 
so far considered only for massless particles. 

The first motivation of this paper is to construct tree-level recursion relations which 
naturally incorporate massive particles. The recursion relations follow from general quan- 
tum field theory arguments, and are valid for any quantum field theory (and also for 
gravity). The advantages of the recursion relations and the MHV rules approaches over 
standard Feynman-diagram calculations are obvious. Just as their massless parents [3, 4], 
the recursion relations for massive particles do not suffer from the factorial growth of the 
number of contributing Feynman diagrams. Compared to MHV rules, the recursion rela- 
tions give more compact results. They can be rigorously derived for a general quantum 
field theory and they incorporate massive particles in a natural and universal way. Essen- 
tially, with the general recursion relations in place, one may be able to avoid in future all 
Feynman-diagrams calculations of nontrivial tree amplitudes. 

In this paper we will use massive recursion relations to derive compact amplitudes for 
gluons coupled to massive coloured scalars. In the forthcoming companion paper [24] we 
will use the recursion relations for amplitudes with massive particles with spin. 

At one loop, the MHV rules have been successfully applied to supersymmetric am- 
plitudes [25-31]. New classes of amplitudes have been derived with the new methods in 
Refs. [5, 32-42]. All of these results give complete amplitudes only in supersymmetric 
theories which are cut-constructible in 4 dimensions [6, 7]. In non-super symmetric gauge 
theories, the new methods apply only to the 4D cut-constructible parts of the amplitudes. 

The second motivation of this paper is to assemble together the pieces necessary for 
the complete calculation of one-loop amplitudes in non-supersymmetric gauge theories. 

One-loop multi-gluon amplitudes in the 'all-plus', and in the 'one-minus' helicity config- 
urations, (+, +,+,...,+) and (— , +,+,..., +), are known [38, 43, 44]. However, the full 
one-loop amplitudes even in the simplest non-supersymmetric gauge theory, pure Yang- 



- 1 - 



Mills, are still known only for n < 5 external gluons, with the five-gluon amplitude cal- 
culated in 1993 in Ref. [45]. The six-gluon amplitude has not yet been calculated for all 
helicity configurations. One- loop amplitudes with n external gluons in pure Yang-Mills are 
conveniently decomposed as 

^gluon J^ = ^ ^^chiral M=\ _j_ ^scalar Q 

where the first and the second terms are the contributions of the M = 4 supermulti- 
plet, and the chiral jV = 1 supermultiplet running in the loop. These contributions arise 
from supersymmetric theories, they are cut-constructible and largely known (at present 
all colour-ordered sub-amplitudes are known for A%= 4 for n < 7, [32, 33, 37], and the 
complete A^ ral N=l is known for n < 6, [6, 7, 30, 31, 34, 35, 39, 42]). 



The last term in ( |1.1|) corresponds to an n-gluon amplitude with a complex scalar 
propagating in the loop. This is a non-supersymmetric one- loop amplitude and it is not 
cut-constructible [7]. What this means is that the amplitude .A^ calar cannot be fully recon- 
structed from its imaginary part evaluated in 4 dimensions. In other words, the knowledge 
of the cuts of the amplitude is insufficient to recover the full answer. The answer contains 
purely rational terms, these do not have cuts in 4 dimensions. 

To regulate infrared and ultraviolet divergent integrations over loop momenta, loop 
amplitudes are commonly evaluated in D = 4 — 2s dimensions. In D dimensions (with 
non- integer D), the amplitudes contain only cut-constructible contributions [46, 47]. This 
is because all terms in the result are proportional to independent powers of the kinematic 
invariants, which give rise to terms like (— s)~ £ that necessarily contain logarithms [46]. 
Hence, in D-dimensions all amplitudes are completely determined from their cuts and 
are cut-constructible. The price to pay is that we now need to know on-shell tree-level 
amplitudes in D-dimensions and work to higher order in e. 

In order to determine full amplitude .A^ calar we cut the propagators in the loop and 
put them on-shell in D-dimensions. On both sides of the cut we are left with tree-level 
amplitudes of the form 

<A%%2(</>\ 1 > 91,92, 9m, </>h) ■ ( L2 ) 
These are the colour-ordered subamplitudes with two adjacent scalars and m gluons of 
arbitrary helicities. The external gluons have four-dimensional momenta pi, ■ ■ ■ ,p m - The 
scalar momenta l\ and I2 are D-dimensional, they are the loop momenta which we cut. 

D-dimensional massless momenta of scalar fields can be thought of as the 4-dimensional 
on-shell momenta of particles with mass fj,, such that l\ = l\ = fi 2 . Here l\ and I2 are 4- 
dimensional and the u 2 term arises from the extra D — 4 = — 2e dimensions [46, 47]. 

Hence, in order to calculate 6-gluon one-loop amplitudes we need to know tree-level 



4, 5, 6-point amplitudes (1^2) with m = 2,3,4 gluons for all independent helicity arrange- 



ments. In Ref. [47], the amplitudes (1^) were calculated for up to 4 gluons with the same 
helicity. Using the massive recursion relation we will reproduce these results as well as 
derive the remaining amplitudes (|1.2|) with positive and negative helicity gluons. 
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Notation 

We will be using the spinor helicity formalism [48] . A vector is represented as a bispinor 

Paa = °aaP» . ( L3 ) 

where are the chiral gamma matrices. The norm of p is PuP^ = det(p aa ). Hence, a 
vector is null if and only if it can be written as a product of two spinors 

p 2 = p ah = \ a \a ■ (1.4) 

Here, A a and A a , are commuting spinors of negative and positive helicity respectively. For 
real momenta in Minkowski signature A is a complex conjugate of A, hence A is convention- 
ally called the "holomorphic" and A is the "anti-holomorphic" spinor. When discussing 
the recursion relations, it will be necessary to consider complex momenta. For complex 
momenta, A a and A a are independent complex variables. 

We will often use a short-hand notation: 

K = \i) a K = W- (1-5) 
The invariant spinor products are defined as 

(i 3) = = VA, a , [i j] = <* + |r> = \% a ■ (1-6) 

Here spinor indices are raised and lowered with e-symbols. The scalar product of a null 
vector pi with a vector q is 

Pi -Q = ~\ • (1-7) 

Throughout the paper we use the sign conventions 1 of [1-3] and define 

Sij...k = (pi+Pj + ...+ Pkf (1-8) 

(i\M = (i\ a Pa h \j]a = -KPaa~^ (1-9) 
(i\PrPs\j) = (i\ a Pra h P J \j) b = -KPraa X jb (1-10) 
\i\PrPs\j] = [i\>ra a Psa' \j] b = ~KPrab pf \ b (1.11) 

For massless momenta equations ( |l.8|) -( |i~TT| ) become 

s ij = (ifiiij] , (i\Pr\j] = (ir)[rj] , (1.12) 
(i\PrPs\j) = (ir}[rs}(sj) , [i\p r Ps\j] = [ir](rs)[sj] , (1.13) 

and so on. 



1 Note that the opposite sign convention from ( ]l.6| ) is used in QCD literature for [i j}. 
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2. Recursion relations with massive particles 

Consider a tree level scattering amplitude of n incoming particles, some of which might be 
massive 

A(pi,P2, ■ ■ ■ ,Pn) , Pi = m- . (2.1) 

Single out two particles, for special treatment. These particles can be either massive 
or massless. For given Pi,pj pick a null vector 77 = that is orthogonal to both pi and 

Pi 

V ■ Pi = f] ■ pj = rj 2 = • (2.2) 

For generic pi and pj, there are exactly two such 77 up to scaling. To see this, consider 
the plane spanned by p- L and pj. Geometrically, the first two conditions in eq. ( |2.2| ) mean 
that 77 lies in the plane orthogonal to the plane spanned by Pi,Pj. The last condition sets 
rj to be in the intersection of this plane with the lightcone. A generic plane going through 
the origin intersects the complex lightcone at two complex rays. These rays define the two 
solutions for 77 up to a scaling by a complex number. We now will construct these solutions. 

Solution for the shift momentum 

Let us find now explicit solutions of eqs. ( |2.2| ) for complex momenta. We will discuss 
in turn the case when both i,j are massless, when one of i,j is massive and at last the case 
when both i,j are massive. 

If the marked momenta are null, pi = AjAj, pj = XjXj, then the condition that rj is 
orthogonal to pi gives 2rj ■ pi = {rj, i)[r), i] = 0, so either A^ = Aj or A,, = Aj,. Similarly, 
vanishing of the Lorentz invariant product of r\ and pj implies A^ = Xj or A^ = Xj. 
Combining these two conditions, we find two solutions 

V = XjXi , rj' = XiXj . (2.3) 

Now consider the case where the particle i is massless and the particle j is massive. 
The condition that momentum r\ is orthogonal to Pi gives — Xi or — A^. If A^ — Aj, 
the orthogonality to pj reads 

2 V -p j = Xtp ja& X^= , (2.4) 

hence A° = A» a p^ a = (XiPj) a . So the two possible null vectors orthogonal to both pi and 
Pj are 

rT = XUXiPjf , V' ad = (PjX^Xf . (2.5) 

The case when i is massive and j is massless is treated analogously. 

The last case to consider is when both particles i and j are massive. Here, neither 
Pi nor pj is a product of two spinors so the expression for 77 is not as simple. We use the 
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condition 2r\ ■ pi = A° A^ pi aa = to express X va = pi aa A^ . Putting this into the second 
orthogonality condition gives a quadratic equation for A^ 

X a v X b v p taaPjbb e hb = -(v\Pi- Pj \v) = • (2.6) 

This equation has two solutions, A^, which we can find, for example, by setting X® = (1, x) 
and solving the quadratic equation for x. The analogous condition for the positive helicity 
spinor A^ 

K, ^ Piaa P jbb e ab = -HPi-pM = (2.7) 

has also two solutions. Altogether, up to scaling, there are two null vectors r\ = A^A,, 
that are orthogonal to Pi,Pj- We do not know of a convenient Lorentz invariant solution to 



eq. (2J3) and eq. (2/7). This makes the case where both marked particles are massive less 



tractable than the two simpler cases where at least one of the marked particles is lightlike. 
2.1 Derivation of the recursion relations 



To construct massive recursion relation for a tree- level n-particle amplitude A{p\ , P 2, ■ ■ ■ , p 
we first mark particles i and j for special treatment and pick one of the two null vectors 



r\ satisfying the conditions in eqs. (2.2). Following [4], consider the auxiliary function of 
one complex variable 

A{z) = A( P i(z), . . . ,Pi(z), . . . , P j(z), . . . , Pn (z)) , (2.8) 

where Pk{z) = Pk for k / i, j , and 

Pi(z) = pi + Z7] , pj(z) = pj - zrj . (2.9) 

Since r\ is null and orthogonal to Pi,Pj, the shifted momenta are on-shell 

Pi (z) 2 =p*, p j {zf=p 2 j . (2.10) 

Equations ( fEs| ) imply that Pi(z) +pj(z) = pi +Pj, so A(z) obeys momentum conservation. 
Hence, it is an on-shell scattering amplitude of particles with complex momenta and can 
be computed from the usual Feynman rules. 

Clearly, the momenta of the external particles are linear functions of z. Notice that 
the spinors of massless external particles are linear functions of z as well. In the case 



where both marked particles are massless, there are two possible r/'s given by eq. (2^). For 
r\ = XjXi, the shift (|2.9| ) is accomplished by 

Xi(z) = Xi + zXj , Xj(z) = Xj - zXi . (2-11) 

The second solution for the shift vector, r( = XiXj, gives 

Aj(z) = Xi + zXj , Xj(z) = Xj — zXi . (2-12) 
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Consider now the case when one of the particles, say particle i, is massless and the 
other particle j is massive. Then eq. ( |2.5| ) gives rj aa = X^(XiPj) a . The shift of marked 
momenta fl2.9|) is accomplished by 

(2.13) 



A? (z) = Af + *(A^)« , pf (*) = pf - zXUXiPjf ■ 
For -rf a = (pj\i) a \f there are analogous expressions 

A? (z) = X? + zipfXtf , pf(z) = pf - zipfXifXl . 



(2.14) 



It follows that A(z) is a rational function of z because at tree level, the scattering 
amplitude is a rational function of the spinors of massless external particles and of the 
momenta of massive external particles. 

s s + 1 




r - 1 



Figure 1: Diagrammatic representation of the recursion relation. Arrows label the momentum 
flow. 

At tree-level the rational function A(z) can only have simple poles in z coming from 
internal propagators 1/P(z) 2 . Each propagator divides the external particles into two 
groups, the particles to the 'left' and to the 'right' of the propagator as illustrated in figure 
|l[ Hence, the momentum P{z) of a propagator is the sum of the momenta of the external 
particles to the left of the propagator 

P = Pr + ...+ Pl + ...+ Ps . (2.15) 

Momentum P(z) depends on z only if the particles i and j are on opposite sides of the 
propagator. We choose the particle i to be on the left of the propagator and the particle j 
to be on the right, as in figure |l[ Then 

P{z)=P + zr), (2.16) 

and the propagator is, 

(2.17) 



P{z) 2 — m 2 P 2 — m 2 + 2zP ■ rj ' 



where m is the mass of the internal particle. The propagator (2.17) has a simple pole at 

P 2 -m 2 , 

(2.18) 



2P-7] 
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For generic external momenta, all internal momenta are different, hence, the locations of 
all poles are different. It follows that the tree-level amplitude A(z) has only simple poles 
as a function of z. 

To find the recursion relations, we use the familiar theorem from complex analysis that 
the sum of residues of a rational function on a Riemann sphere is zero. Applying this to 
A{z)/z we express -4(0) as a sum over residues 

,4(0) = Res (*&) =-£Res(^( - Res ( *$) , (2.19) 

V z J 2 =0 a V Z J z=z a V z J Z =oo 

where the sum is over all finite poles z a of the amplitude A(z). These come from the prop- 
agators 1/P 2 (z) that separate the particles i and j. The residues at finite z are determined 
by the factorization of the scattering amplitude when the Feynman propag ator ( [DTD goes 
on-shell 

Res (42) = - ^(za)A R (z a ) 

Here, Al and Ar are the tree-level amplitudes of the particles to the left and to the right 
of the propagator and z a is given by eq. ( 2.18 ). 

Hence, any tree-level scattering scattering amplitude A = -4(0) can be written in the 
form 

A = +r, AL( % )AR i Za) -Res , (2.21) 

„ Pi -mi V z Joe 

where the sum is over all channels a such that the particles i and j are on different sides 
of the channel, and z a is given by eq. ( 2.18| ). The relation (2.21) is useful for computing 
scattering amplitudes only if there is an efficient way to determine the boundary contribu- 
tion Res (A(z)/z) 00 . The most favourable scenario is when this contribution vanishes. This 
happens if and only if A(z) vanishes at infinity, 

Res (A^) =0 A(z) -► 0, for z -> oo , (2.22) 



in which case, there is a simple recursion relation: 

A = Y,A L (z a ) l _ 2 A R {z a ) (2.23) 

that expresses A in terms of lower-point on-shell scattering amplitudes Al(zoi) an d Ar{zq). 
The summation in (2.23) runs over all partitions of particles between Al{ z o) and AR(z a ), 
such that pi is on the left, and pj is on the right of P a , and also over all helicities h of the 
intermediate state P a . 

The above considerations apply to the case with massive or massless marked parti- 
cles. However, for calculations carried out in this paper it will be sufficient to take both 
marked particles to be massless. In this case, the necessary conditions for the vanishing of 
the boundary contribution ( |2.22| ) put constraints on the possible helicities of the marked 
particles i and j. We discuss these conditions in section |4| 
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2.2 Recursion relations: summary 



We will use the recursion relations to calculate tree-level scattering amplitudes in Yang- 
Mills theory coupled to matter fields. The matter fields may be massive or massless and 
transform in a generic representation of the gauge group. We consider the colour-ordered 
partial amplitudes A = A(p\, ■ ■ ■ ,p n ), i n which the coloured particles come in a definite 
cyclic order 1,2, ... ,n. These amplitudes are obtained by stripping away the colour factors 
from the full amplitude, hence, they depend on the kinematic variables, momenta and 
helicities, pk and hk only. 

In the remainder of the paper we will take both marked particles to be massless. We 
shift two massless momenta pi = and pj = \j)\j] of the marked particles by r/ = \j)\i], 
so the shifted momenta are 

Pi =Pi + z\j)\i] , (2.24) 
Pj = Pj - z\j)\i] , (2.25) 
P = P + z\j)\i], (2.26) 

where P = p r + . . . + pi + . . . + p s is the momentum of the intermediate particle. For the 
particles i,j this is equivalent to shifting the spinors 

@ = l»1 > (?) = N> + Ai) - (2-27) 
l?) = |i>, \3} = \j]-z\i] ■ (2-28) 

The recursion relation ( |2.23| ) written more explicitly is (c.f. figure [l]) 

An(pi, . . . ,p n ) = V V A L (p r ,...,Pi,...,p s ,-P h ) — 2 

partitions h p 

xA R (p- h ,p s+ i, . . . ,Pj, . . . ,pr-i) , (2.29) 



where summation is over all partitions of n external particles between Al and Ar, such 
that pi is on the left, and pj is on the right, and also over the helicities, h, of the intermediate 



state, z can be found from the on-shell condition P 2 = m 



P 2 — m 2 P 2 — m 2 

*=-W=l^#' (2 ' 30) 

When the intermediate state P is a massless particle (e.g. a gluon), we can simplify 
the spinor products involving \P) and |P] as in Ref. [3]: 

, ~ (k\P\i] (k\P\i\ 

{k p) = iMi = mi , (2.3D 

[h]^,^, (2.32) 
(j P) w 
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where lo and u enter the amplitude always in the combination ujuj = (j\P\i\. 

For practical computations it is essential that A{z) vanishes for large z, so that the 
recursion relations do not have a boundary contribution at infinity. As discussed in section 
|4], this puts a constraint on the helicities of the particles i and j. For our choice of the shift 
momentum, r/ = \ j)\i], the helicities of the marked particles can take the values, 

V=\j)\i] : (huh,) = (+,-) , (+,+) , (-,-) (2.33) 

but not (hi,hj) = (—,+). Conditions (|2.33| ) are the same for massive and for massless 
amplitudes. 



3. Amplitudes with gluons and massive scalars 

In this section we consider scattering amplitudes of gluons with massive complex scalars. 
These amplitudes are related to amplitudes with massless scalars that have D-dimensional 
momenta. The scalars with D-dimensional momenta Pd can be thought of as massive 
scalars in 4 dimensions. The D-dimensional on-shell condition, Pp 2 = 0, gives the 4- 
dimensional mass-shell equation, P4 2 = /j 2 , where the mass term, /i 2 , arises from the extra 
D — 4 dimensions of momenta. 

We will derive amplitudes with 2 scalars and up to 4 gluons with arbitrary helicity 
configurations. The amplitudes with the same-helicity gluons have been previously derived 
in [47]. 

3.1 Primitive vertices 

The recursion relations construct n-point amplitudes from on-shell m-point amplitudes 
with m < n. The m-point amplitudes are connected to each other with scalar propagators. 
Using the recursion relation n — 3 times gives a representation of the n-point amplitude 
entirely in terms of the 3-point vertices. Hence, 3-point vertices are the building blocks of 
the amplitudes in the recursive approach 2 , they will be called the primitive vertices. 

The recursion relation reduce the task of computing general amplitudes to the compu- 
tation of all 3-point primitive vertices. In this paper we consider amplitudes with massless 
gluons g and massive scalars (j). These can be built from the ggg and (fjgcfi vertices. 



The three-gluon primitive amplitudes can have h or + -| — helicity configurations. 

These are the standard MHV and MHV 3-point on-shell amplitudes 

M9i,92,9t) = , Mgf, 91,9s) = • (3.1) 



2 In particular, this implies that the 4-point vertices in the microscopic Lagrangian are not used in the 
recursive construction of gauge-invariant amplitudes [3]. 



-9- 



The gluon momenta k{ are assumed to be complex which ensures that these amplitudes do 
not vanish on-shell [1, 3]. 



In order to compute scattering amplitudes of gluons and massive scalars, we need to 
determine the (j>g<ft vertices. To obtain these, we start with the off-shell Feynman vertex 
of two scalars of mass /U and momenta li, I2, and a single gluon with momentum k, 

V 3 (lt,k»^)= ^ • (3-2) 

The V2 comes from the normalization conventions used in colour-ordered Feynman rules 
[49], and the + and — indices are labels for a scalar and an anti-scalar. To derive the desired 
on-shell amplitudes, As(lf ,k ,1^), we contract Vs(lf , k 1 *, 1%) with the gluon polarization 
vector, e^(k, q) 



e+(k,qU = V2*d* , e-(k, Q ) = - v^|| , (3.3) 

where q = \q)\q] is an arbitrary reference vector that is not proportional to k. The two 
independent on-shell vertices immediately follow 

A 3 (l+,k + ,l;) = A 3 (l^,k + ,l+) = -M^l , (3.4) 

A 3 (l+,k-,q) = Mll>k-,l+) = "-^p • (3-5) 

We have already noted that the primitive vertices vanish for on-shell real momenta 
in Minkowski space, but are nonzero for on-shell complex momenta. Indeed, the on-shell 
conditions, l\ = l\ = an< i k 2 = 0, together with the momentum conservation imply that 
the momentum of the gluon is orthogonal to the momenta of the scalars k-l\ = k-l% = 0. For 
real massless momentum k in Minkowski space, the spinors k a , ka are complex conjugates 
k* = ±/c a . Similarly a real massive momentum forms a Hermitian matrix (I l)' = lha- 
Hence for real momenta, the conditions lf a k a ka = 0, i = 1,2 imply lf a k a = lf a ka = 
for i = 1,2. It follows that the 3-point vertices (3.5-3~4j) vanish. For example we have 
A(lt,k+,q)ocq a lfk = 0. 

For complex momenta the spinors k a , k a become independent variables. This additional 
freedom allows us to take the momenta of the scalars on-shell while keeping the three- valent 
amplitudes nonzero. 

Finally, we note that the primitive amplitudes are gauge-invariant, even though eqs. ( |3.4| ) 
Q3.5| ) contain explicit ^-dependence. Different choice of the reference vector q amounts to 
a gauge transformation, hence the on-shell amplitudes should not depend on the choice of 
q by virtue of gauge symmetry. It is easy to see this explicitly e.g. for the (4> + g + (f)~) am- 
plitude. The reference spinor q a , a = 1,2 lives in a two dimensional complex vector space. 
The spinors q a and k a are independent due to the condition (qk) 7^ 0, so we take them as 
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a basis of the vector space. Hence a change in the reference spinor can be parameterized 
as 

q ' a = aq a + f3k a . (3.6) 

Changing q, the amplitude becomes 

A'n+ h+ j-\ - a (q\h\k\ + [3 {k\h\k} 

A^,k ,l 2 )- ^ . (3.7) 

Here, = —2k ■ l\ = fi 2 — l 2 = is zero by momentum conservation. The remaining 

a dependence gets cancelled between the numerator and the denominator leaving us with 
the original amplitude. 

It follows that the choice of the reference momenta q% of the gluons does not affect 
the amplitude so in principle we could set them to arbitrary values. In the following 
sections, when using recursion relations to calculate amplitudes with scalars, we will find 
it convenient to set the reference momentum of a marked gluon in a primitive vertex ( |3. 4j ) 
or ( |3.5| ) to be the momentum of the other marked gluon. 

In the following sections we will calculate tree-level amplitudes of the form 

-4(^,31, 02> ... ,g m , <f>i 2 ) ■ (3.8) 

These are the colour-ordered subamplitudes with two massive scalars and m gluons (2 < 
m < 4) of arbitrary helicities. 

When scalars transform in the fundamental representation of the gauge group, the 
'string' of fields in the amplitudes must always start and end with the scalar, precisely 
as in fl3,8|) . Using cyclic symmetry of colour-ordered amplitudes, the scalars in (|3.8| ) can 
be thought of as adjacent. Scalars in the adjoint representation can appear anywhere in 
the string, i.e. they do not have to be adjacent. Such amplitudes can also be calculated 
straightforwardly with our methods. 



We will determine all the independent helicity configurations in ( |3.8| ), all the remaining 
configurations can be obtained from those via the following identities: 

A m+2 {<Pt , gl 1 , g h 2 2 , • • • , gt , ^) = A m+2 , g^ ,g h 2 \..., gt > 4>t 2 ) (3.9) 

= (-ir^ m+2 (^,^,..., 52 h2 , 5l fel ,<A+) (3.10) 
= A* m+2 (^ i ,g^ h \g 2 h ^...,g- hm , ( f>+) (3.11) 



where . . . indicate gluon fields and hi is the helicity of the i th gluon. Equations ( 3.10| ) and 



( 3. lip follow from reflection and parity symmetry of the colour-ordered amplitudes, and 



( pO| ) follows from eqs. 
3.2 4-point amplitudes 



There are two independent helicity amplitudes in this case, the recursion relation (gjp 
gives only one term for each of the amplitudes, as illustrated in figure |2[ 
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Figure 2: Representation of the 4-point amplitude using the recursion relation with 1 and 2 as the 
shifted momenta. 

In the case of two positive helicity gluons, we have: 

Mil, 2 + , q) = A 3 (lf,l + , -P-) 1 2 A 3 (P + ,2 + , l 2 ) , (3.12) 

r n 

where we took the marked particles to be the gluons with momenta k\ and k 2 - We set the 
reference vectors q\ and q 2 of the two gluons equal to the marked momenta gluon on the 
opposite side of the diagram 

q 2 =k 1 = \ki)\kx\ , qi=k 2 = ^2)^2} ■ (3.13) 

We shift the momenta along the vector r\ = |2)|1], so that |1] = |1] , |2) = |2). Whence, 
the amplitude becomes 

(gjjijl] 

Using l\ + l 2 + k% + k 2 = = l\ + l 2 + k% + k 2 , this can be written as 

trft^ii) (12) [12] (Ji • h) - 2(h ■ h)(k 2 ■ h) 



2 (I2) 2 ((h + k,y - fi 2 ) (12)2 ((Zx + k{) 2 _ ^2) 

To get the second expression we used a Fierz identity. The second term in the numerator 
vanishes, k\ ■ l\ = 0. The easiest way to see this is to use momentum conservation in the 3- 
point vertex, l\ + k\ = P, and the on-shell conditions l\ = P 2 = fi 2 , k\ = 0. Alternatively 
this can be shown with the use of the definition fl2~2l ) of k 

= 2(2 | Zl | l3 (SUxIl] = ■ 

This leaves us with the final answer 

Mtt,i*X,H)= (12)( / + [ g,_, v (3-14) 

This agrees with the previously known result computed by Bern, Dixon and Kosower [47]. 

For the amplitude with one positive helicity gluon and one negative helicity gluon the 
recursion relation in figure |2] yields, 

A±(lf,l + ,2~,l 2 ) = A 3 (lf,l + ,-P-) 1 A 3 (P + ,2~,q) . (3.15) 

P l — jJL Z 
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Using the same choice for the marked gluons and the reference vectors ( 3.1 3| ) as before 
gives the result 

<2|/i|l] 2 



A4(l± ,1 ,2 ,l 2 



(12) [12]((Z 1 + £4)2-^) 



(3.16) 



which we checked against a Feynman diagram calculation. 



3.3 5-point amplitudes 

The amplitudes with three gluons and a pair of scalars have three independent helicity 
configurations. As before, we mark the gluons with momenta k\ and &2, and pick their 
reference momenta to be q% = %2, q% = k\. The recursion relation is depicted in figure 
|3[ For the amplitude with all gluons of positive helicities, the recursion relations have 



It ~ 

Z+ 1+ It 1+ l 2 , 1+ 




Figure 3: The decomposition of the 5-point amplitude using the recursion relation with 1 and 2 
as the shifted momenta. 

single non-zero diagram. The diagram with gluon exchange vanishes as the choice of shift 
vector tj implies vanishing of the A{2 + , 3 + ,p ) MHV amplitude. The amplitude follows 
immediately: 

A (1+ 1+ 2+ 3+ D - ^ [3|(l + 2)/i|l] , 17 , 

5il ' ' ' , 2 '~ {{h + kx? -tf) (12) (23)({h + ■ [6An 

This is in agreement with the result in [47]. 

For the case where one of the gluons has negative helicity we have two independent 
helicity configurations, each of which has two non-zero contributions: 

(3|/ 2 (l + 2)/ 1 |l] 2 



•^5(^1 j 1 i 2 ,3 , Z 2 



((h + k,) 2 - M 2) (12) (23)((Z 2 + k 3 ) 2 - (j,*) [3|(1 + 2)h\l] 
f [12] 3 



s 12 3 [23] [31(1 + 2)^11] ' ( ' ! ' LS) 



((Zi + A;i) 2 - /i 2 ) (12) (23)((Z 2 + Zc 3 ) 2 - fi 2 ) [3|(1 + 2)fc|l] 

+ ^ t 13 J 4 (3 19) 

^ s 123 [12] [23] [3|(l + 2)i!|l] • { - ) 
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These results are new. Our results ( |3.17 )-( 3.19| ) numerically agree with the much lengthier 
expressions which we obtained by a direct calculation of the 25 Feynman diagrams. 



3.4 6-point amplitudes 



We mark gluon momenta 1 and 2, and write down the recursion relation for the 6-point 
amplitudes with 4 gluons in figure ^. 




V 



+ 



i; 



\ / 
\ / 



It 



1+ 



Figure 4: The decomposition of the 6 point amplitude using the recursion relation with 1 and 2 
as the shifted momenta. 



In the case of all gluons of the same helicity, only the first diagram contributes. We 



find, 



A(^,l + ,2+3+,4+ Z" 



fi 2 [4|Z 2 (3 + 4)(l + 2)fi|l] 
Q1Q2Q3 (12) (23) (34) ' 



where Q 1 = ((h + h) 2 - fi 2 ), Q 2 = {{h + h + k 2 ) 2 - fi 2 ) and Q 3 = ((Z 2 + k 4 ) 2 - /U 2 ) 



Eq,3.20 is a slightly shorter form of the result given in [47]. 



Now we compute the remaining independent 6-point amplitudes. There are two am- 
plitudes with one negative helicity gluon: 



A(/r,l + ,2 + ,3+,4-,Z 2 -) = 

(Q2 (4JM1 + 2 + 3)fi|l] - fi 2 ( <4|Z 2 32|1])) 2 
Q1Q2Q3 (12) (23) (34) [4|/ 2 (3 + 4)(l + 2)/ 1 |l] 

, ^ [3|(l + 2)/ 1 |lf 

+ Q 1 (12) [34] (2|(3 + 4)(h + l 2 )h\l] [4|/ 2 (3 + 4)(1 + 2)^1] 

li 2 (4|2 + 3|1] 3 

51234^234 (23) (34) <2|(3 + 4)(Zi + J 2 )ii|l] 
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A(^,l + ,2 + ,3-,4+,/ 2 -) = 

(Q 2 (31/illl-zi 2 [32] [21]) 2 (3|/ 2 |4] 2 



+ 
+ 



Q1Q2Q3 (12) (23) (34) [4|/ 2 (3 + 4)(l + 2)Z 1 |l] 

[4|(1 + 2)Z 1 |1] 4 



Qi (12) [34] [3|(1 + 2)ii|l] (2|(3 + 4)(h + h)h\l] [4|Z 2 (3 + 4)(1 + 2)Zi|l] 
/i 2 (3|2 + 4|1] 4 



si 23 4S 23 4 <4|2 + 3|1] (23) (34) (2|(3 + 4){h + l 2 )h\l] 
[4|(1 + 2 + 3)Z 1 |l] [Iff 
+ Sl23 (4|2 + 3|l]Q 3 [23][3|(l + 2)/ 1 |l] 1 ' j 

These amplitudes agree with the massless MHV-type amplitudes as /i 2 — > 0. There are 
three independent helicity amplitudes with two negative helicity gluons: 

A(^,l + ,2 + ,3-,4-,/ 2 ") = 

(Q2 (4|Z 2 (3 + 4)Zi|l] - M 2 (4|/ 2 (3 + 4)2|1]) 2 



+ 
+ 
+ 



Q1Q2Q3 (12) [34] (4|/ 2 (3 + 4)|2) [3|(1 + 2)Zi|l] 

f (2|/i|l] 2 (34) 3 

Q 1 (12) (23) <2|(3 + 4)(Z 1 + Z 2 )Z 1 |1] (4|/ 2 (3 + 4)|2) 

[l|Mi|l] 2 (34) 3 

S1234S234 <4|2 + 3|1] (23) (2|(3 + 4)(/ 1 + Z 2 )/ 1 |l] 

(4|^(l + 2 + 3)^|l] 2 [12] 3 

Q3S123 <4|2 + 3|1] [4|(l + 2 + 3)/i|l] [23] [3|(1 + 2)Zi|l] 
[}} [12] 3 

+ ^1234 [23] [34] [4|(l + 2 + 3)Zi|l] (3-22) 



AGi + ,l + ,2-,3 + ,4-,/ 2 -) = 

(2|/i|l] 2 (2|/ 1 + 1|3] 2 (4|l 2 |3] 2 

Q1Q2Q3 (12) [34] (4|/ 2 (3 + 4)|2) [3|(1 + 2)/i|l] 
f (24) 4 (21/xll] 2 



+ 



Qi (12) (23) (34) (2|(3 + 4)(h + l 2 )h\l] (4|/ 2 (3 + 4)|2) 

[IK2/1II] 2 (24) 4 

S1234S234 (4|2 + 3|1] (23) (34) (2|(3 + 4){h + l 2 )h\l] 

(4|/ 2 (l + 2 + 3)/ 1 |l] 2 [13] 4 

g 3 si23 [12] [23] [4|(l + 2 + 3)Zi|l] [3|(1 + 2)/i|l] <4|2 + 3|1] 

+ (3 23) 

+ Sl234 [12][23][34][4|(l + 2 + 3)/ 1 |l] 1 ' ) 
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A 6 (lf,l + ,2-,3-A + ,l 2 ) = 

(2|/ 1 |1] 2 (Q 3 [4|f 1 + l|2)-/i 2 [43] (32)) 2 
Q1Q2Q3 (12) [34] [3|(l + 2)Zi|l] (4|/ 2 (3 + 4)|2) 

^ 2 (23) 3 (2\h\lf 

Qi (12) (34) (4|/ 2 (3 + 4)|2) (2|(3 + 4){h + l 2 )h\l] 

[l|Mi|l] 2 (23) 3 

S1234S234 (34) (2|(3 + 4)(/ 1 + Z 2 )Z 1 |l] (4|2 + 3|1] 

[l|(2 + 3)Zi|l] 2 [l|(2 + 3)/ 2 |4] 2 

S123Q3 [12] [23] [3|(l + 2)Zi|l] [4|(l + 2 + 3)/i|l] <4|2 + 3|1] 



si234 [12] [23] [34] [4|(l + 2 + 3)/i|l] v " ' 

All of the above six-point amplitudes have been checked numerically against an independent 
calculation of the 220 Feynman diagrams. 



4. Vanishing of A(z) at Infinity 



The recursion relations are valid as long as A(z) — ► for z — > 00. This is the case only for 
some choices of marked particles. Indeed, for the marked gluons i,j with 

\i) = \i) + z\j) , \j] = \j]-z\j] , (4.1) 

it was observed in [3] that A(z) — > for the helicity assignments of gluons (hi,hj) = 
(+,+),(+,—),(—,—). This has been extended to fermions in [19, 20]. A direct diagram- 
matic proof for gluon amplitudes in the (+, — ) case was given in [4]. Here, we present an 
argument valid for all three helicity assignments of marked gluons putting no restrictions 
on the remaining partons 3 . We also consider the case where one or both of the marked 
particles are massless fermions. 



4.1 The (+,-) case 



We start with the (+, — ) case for which we will show that all Feynman diagrams con- 
tributing to A(z) vanish at infinity. In a given Feynman diagram, the z-dependence flows 
along a unique path of Feynman propagators and vertices. The Feynman diagrams with 
most dangerous z-dependence are those in which all vertices along this path are trivalent. 
These vertices are linear in momentum so they each contribute a factor of z. Each of the 
propagators gives a factor of 1/z. A path made of r propagators has r + 1 vertices. Hence, 
the propagators and vertices go at infinity as z r+1 / z r = z. For fermions, the propagator 
(f+m)/{p 2 — m 2 ) goes like a constant at infinity, but so do the vertices containing fermions. 

3 For the (+,+) case, our proof works for the amplitudes that have two extra negative helicity gluons. 
However, see the end of this section for a heuristic argument, which removes this assumption. 
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Each of the fermion propagators increases the above estimate by a factor of z and each of 
the fermion vertices decreases the estimate by z. For / fermion propagators, there are at 
least / + 1 fermion vertices along the path of the z-dependence, hence the above bound on 
the Feynman diagram gets improved by at least a factor of 1/z. 

The remaining pieces of the Feynman diagram that depend on z are the polarization 
vectors of the marked gluons. Consider marked gluons of opposite helicity, so that we can 
have either gf, gj , or g~ , gj. The spinors \{(z) and Xj(z) are linear in z, eq. ( 2.11] ), while 
\i(z), Xj(z) are independent of z. It follows from eq. ( |3.3[ ) that for the helicity configuration 



9i >9j ! the polarization vectors of both gluons go as 1/z at infinity. Altogether Feynman 
diagrams vanish at infinity as 0(z/z 2 ) = 0(l/z). 

In the opposite case, g~ , g^, gluon polarization vectors grow as z at infinity and 
individual Feynman diagrams go as C(z 3 ) leading to nontrivial boundary contributions to 
the recursion relations for A{z). 

If one of the marked particles is a gluon and the other a massless fermion, we can 
proceed analogously. Now the propagators and vertices contribute a factor of at worst z° 
because of the fermions. For helicity assignments (hi,hj) = (+1, —1/2), (+1/2, —1), the 
gluon polarization vector goes like 1/z and the fermion polarization spinors u~(pj) = A" 
and u + (pi) = A? are independent of z. Hence the shifted amplitude vanishes at infinity for 
the assignment (hi,hj) = (+,—). 



4.2 The (+,+) and (-,-) cases 



The (hi, hj) = (+, +) case cannot be treated by counting powers of z in individual Feynman 
diagrams and requires a more elaborate argument along the lines of [22]. The argument 
works for amplitudes that, besides gluons gf-,9~j have two extra negative helicity gluons 
g^,g^, hence in particular it works for all QCD gluon amplitudes since at tree level, the 
gluon amplitudes always have at least two gluons of both helicities. 

Consider the function A(\ m (z), \ m (z)) constructed from the scattering amplitude A 
by shifting the spinors 

|») = \i) +z\k) +z\l) 
\k] = \k] - z\i] 

\l] = \l] - z\i] (4.2) 

while keeping \k) and |/) unshifted. The shifts ( f4.2j ) maintain momentum conservation, 
so A(z) can be computed from usual Feynman diagrams. 

Note that the function A(z) vanishes at infinity. As in the previous case, we prove 
this by studying the most dangerous Feynman diagrams. The z-dependence flows along 
a 'three-legged path' of Feynman propagators and vertices, with each leg ending at one 
of the external gluons i, k, I This path is illustrated in figure ||. For a path consisting of 
r propagators, there are r + 1 vertices. For large z, the propagators and vertices give a 
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z r+1 1 z r = z contribution to the z-dependence of A(z). The polarization vectors of each of 
the gluons give a factor of 1/z so the function A(z) vanishes for large z as 0(1/ z 2 ). 

Hence, following the discussion of recursion relations in section |2|, we have 

A = ^A L {z a , -Pa{z a ) h ) 2 - 2 A R (z a , P a , (4.3) 

a,h ex. m a 

where the summation is over partitions of the particles into two sets L,R such that i £ L 
and at least one of j, k is in R, and also over helicities of the intermediate state P a . 



The amplitudes Al,Ar are evaluated at 



p2 _ 2 

(4.4) 



(k\P a \i] + (l\P a \i] 



Similar recursion relations with three or more marked particles have been recently consid- 
ered in [21, 22, 38]. 

„•+ 

k- 




Figure 5: In a Feynman diagram with marked particles i,k,l, the z dependence flows along a 
unique three-legged path of propagators and vertices. 

We will use this expression for A in an induction proof of the recursion relations with 
(hi,hj) = (+,+). We consider the function A(y) with 

\i) = \i) + y\j) \j] = \j]-y\i] (4-5) 



and estimate the y dependence of A(y) with the help of eq. ( |4.3j ). For a partition a such 
that i,j are both in L, neither P a nor z a depend on y. So the entire y dependence comes 
from Al- But Al has fewer particles than A so by induction hypothesis it vanishes as y 
goes to infinity. 

The other case with i G L, j £ R is more subtle. We have P a (y) = P a + yXjXi so Z\ 
depends linearly on y 

_ P* - y {j \P a \i} - ml 

(k\p Q \i}+ (l\p a \i] ■ 

We represent the on-shell amplitudes Al,Ar as a sum of Feynman diagrams. The y 
dependence flows along a four-legged tree with each leg of the tree ending at one of the 
gluons i,j,k,l. A tree composed of r propagators (one of which is the propagator connecting 



- 18 - 



Al and Ar), contains r+1 vertices. The most dangerous contribution comes from Feynman 
diagrams in which the r + 1 vertices are three-point vertices that each give a factor of y. 
Hence, the propagators give a factor of l/y r and the vertices a factor of y r+1 . 

The polarization vectors of the gluons i, k, I each give a factor of 1/y and that of the 
gluon j a factor of y, so altogether we have 1/y 2 coming from the polarization vectors. 
There are two more 'external' particles with y dependent momentum. These come from 
the ends of the internal line connecting Al and Ar. They contribute J2h e>le ~ h which, in 
a convenient gauge, goes like a constant at infinity. For example for a gluon we have 

£ ^ = 9r ~ ^Vffr , (4.7) 
h=± q 

where P = P a {y) and q is the reference momentum of the internal gluon. For large y, this 
is independent of y. Hence, A(y) vanishes as 0(y r+1 / (y r y 2 )) = 0(l/y) for large y. 

In this proof, we have assumed that the amplitude has besides the two marked positive 
helicity gluons i,j also two negative helicity gluons k,l. However we believe that this 
restriction is not necessary. Heuristically, one can argue as follows. Consider any amplitude 
A without imposing the restriction that it contains two negative helicity gluons. We can 
construct an 'auxiliary' amplitude A' which, besides external particles present in A also 
has two extra negative helicity gluons k, I. In the previous argument, we demonstrated that 
A'{y) with marked gluons i,j vanishes for large y as 0{\/y). One can imagine taking soft 
gluon limit pi,Pk — > to recover the original amplitude. Taking the soft limit gives back a 
universal factor, times an amplitude with the gluon removed. The universal factor depends 
on the momentum of the soft gluon and of the possible adjacent particles. Hence, as long 
as the 'auxiliary' gluons k, I are not adjacent to the marked gluons i, j, their soft factors are 
y independent. So the y asymptotic of the amplitude A with gluons k, I removed should 
be the same as that of A'. 

5. Summary 

In this paper we have constructed tree-level recursion relations valid in any quantum field 
theory which naturally incorporate massive particles. As an application, we focus on scalar 
particles. However, there is no obstacle in applying our general method to massive particles 
with spin. With the general recursion relations in place, one may be able to avoid in future 
all Feynman-diagrams calculations of nontrivial tree amplitudes. 

As a first application of these recursion relations, we have derived expressions for 
scattering amplitudes involving a pair of massive scalars and up to four gluons. Explicit 
results are given in Section 3. These compact analytic expressions for amplitudes may be 
useful in deriving the one-loop six gluon amplitudes in QCD using the unitarity method of 
Bern et al [46, 47]. 
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